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The equations governing weakly nonlinear modulations of A'^- dimensional lattices are consid- 
ered using a quasi-discrete multiple-scale approach. It is found that the evolution of a short wave 
packet for a lattice system with cubic and quartic interatomic potentials is governed by generalized 
Davey-Stewartson (GDS) equations, which include mean motion induced by the oscillatory wave 
packet through cubic interatomic interaction. The GDS equations derived here are more general 
' than those known in the theory of water waves because of the anisotropy inherent in lattices. Gener- 

alized Kadomtsev-Petviashvili equations describing the evolution of long wavelength acoustic modes 
in two and three dimensional lattices are also presented. Then the modulational instability of a iV- 
dimensional Stokes lattice wave is discussed based on the A^-dimensional GDS equations obtained. 
Finally, the one- and two-soliton solutions of two-dimensional GDS equations are provided by means 
of Hirota's bilinear transformation method. 
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I. INTRODUCTION 



Since the pioneering work of Fermi, Pasta and Ulam |T] on the nonlinear dynamics in lattices, the understanding 
of the dynamical localization in ordered, spatially extended discrete systems has experienced considerable progress. In 
particular, one-dimensional(lD) lattice solitons, which are localized nonlinear excitations due to the balance between 
nonlinearity and dispersion, are shown to exist ||^. Similar to the cases in fluid physics and nonlinear optics, most 
of the analytical approaches on lattice solitons are based on weakly nonlinear theory. The basic idea of the weakly 
ij-^ _ nonlinear theory is that linearized lattice equations are assumed to provide a satisfactory first approximation for those 
\^ . finite-amplitude disturbances which are, in some sense, sufficiently small. Successive approximations may then be 
' developed by an asymptotic expansion in ascending powers of a characteristic wave amplitude. The weakly nonlinear 
, theory has been shown to be very successful in revealing many important physical processes, e. g. resonant wave-wave 
' interactions, modulational instability, the formation of solitons, etc, in a clear-cut way. A very useful method for the 
asymptotic expansion is the method of multiple-scales, which in the case of lattices reduces the system to a set of 
partial differential equations for the envelope (or amplitude) while the original system is a set of differential-difference 
equations, and usulally cannot be solved exactly. There are two basic advantages of the multiple scale expansion: (i) 
. it contains a unique explicit small parameter, and hence is controlable, and (ii) it allows obtaining solutions in an 
^ ' explicit form. It is well known that, for a ID lattice wave with a large spatial extension, the envelope of the lattice 
Q , wave is governed by the nonlinear Schrodinger (NLS) equation for a short wavelength packet [|| and the Korteweg-de 
O Vries (KdV) equation for a long wavelength acoustic mode Q. 

_ ^ In recent years, much attention has been paid to coherent structures in multi-dimensional lattices (see e.g. Ref. 

. [^). In particular we mention a generalization of the KdV equation in a 2D lattice with only a cubic interatomic 
' potential, i. e. the Kadomtsev-Petviashivili (KP) equation, derived for a lattice wave traveling in a given direction |^ 
5^ ] and coupled 2D NLS equations describing quadratic solitons due to the second harmonic generation in a 2D lattice 
of the two-component dipoles 0. However, to the best of our knowledge, up to now 2D and 3D generalization of the 
NLS equation with a mean motion induced by oscillatory wave packets in lattice systems (i.e. due to long wavelength 
acoustic mode) has not been developed. Meantime such motion introduces dramatical changes in the lattice dynamics. 

In the present paper, using a quasi-discrete multiple-scale approach PJ^-p^, we derive generalized Davey- 
Stewartson (GDS) equations in multidimensional lattices with cubic and quartic interatomic potentials. Because 
of the anisotropy inherent in lattice systems (i. e. without continuous translation and rotation symmetries), in the 
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case of two dimensions the GDS equations presented here are more general than those obtained in water waves , 
which are physicaUy isotropic. We also derive a generalized KP equation governing the evolution of a long wavelength 
acoustic excitation traveling in any direction. 

The organization of the paper is as follows. In Sec. II we formulate the model and deduce the equations for 
slowly varying amplitudes in an A'^D lattice. In Sec. Ill we concentrate on excitations in a 2D lattice. The dynamical 
equations in the long-wavelength limit are presented in Sec. IV. In Sec. V we discuss modulational instability of a 
A^-dimensional Stokes lattice wave on the basis of the GDS equations. Section VI provides some one- and two-soliton 
solutions for the 2D GDS equations based on Hirota's bilinear transformation method and demonstrate the effect of 
anisotropy on the soliton formation. The outcomes are summarized in the last section. 



II. MODEL AND ASYMPTOTIC EXPANSION 

The system we consider is a monoatomic scalar lattice with nearest-neighbor interatomic interactions. The 
equations of motion describing the system are given by 

d 

^ J2j[u{n + a..j) + u{n - aj) - 2u(n)] 

d 

+ H J3j{[u{n + aj) - u{n)f - [u{n - aj) - u{n)f} 

d 

+ J2 JijiH^ + aj) - t^(n)]' + [u{n - aj) - u{n)]^}. (1) 



Here u(n) is the displacement from its equilibrium position of the particle having the mass M and located at 
the site n — ^j^i^jaj, nj being integers, aj being the lattice vectors, and d being the dimension of the lattice, 
Jaj = Kaj/M (a = 2, 3, 4), K2.j, K^j and K^j are harmonic, cubic and quartic nearest-neighbor force constants, 
respectively. Notice that the anisotropy of the lattice is included in the consideration (i.e. in a generic case K^j ^ K^^i 
for i ^ j). We include the cubic potential here since most of realistic interatomic potentials (such as the potentials of 
Born-Mayer- Coulomb, Lennard- Jones, Morse, Toda, etc.) display strong cubic nonhnearity (i. e. Ja,3 7^ 0) In 
the most direct physical applications (namely, to atomic crystals) the dimension d can be either 2 or 3, although more 
formal lattices with d being bigger than 3 are available. In the present section we deal with the last, more general, 
case. 

In order to investigate weakly nonlinear modulation of a lattice wave packet we use the quasi-discrete multiple- 
scale method p|P-p^] to derive the envelope equations describing the development of the modulation of the packet 
along the line of Davey and Stewartson for water waves Namely, we set 

M(n) =^e'^7.,(r,r;0(n,t)) (2) 
i/=i 



with 



en 



vt). 



(n, t) = q • n — uit, 



(3) 



where e is a formal small parameter representing the relative amplitude of the excitation, q is the wave vector: 
q = Ij^j^ being the vectors of the reciprocal lattice: • aj = Sji and uj is the frequency of the respective 

harmonic. The constant vector v as well as the link between u and q, i.e. the dispersion relation, are to be determined 
by solvability conditions. 

There are some comments to be made here. In a generic case the entries of the expansion (|^) depend on the whole 
hierarchy of "slow" variables, i.e. one should consider the set of variables {r,^, t^,} {ly = 1,2, ...) where — e'^(n — Vj^t) 
and — e'^t, which are regarded as independent. In the case one is interested in the effect of quadratic and cubic 
nonlinearity only the scales up to r2 and t2 turn out to be relevant. In the present paper we restrict our consideration 
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to the solutions independent on V2 and this a reason of introducing only the "lowest order" slow variables r = ri and 

Substituting (||) and into Eq.(|l]) and equating the coefficients of the same powers of e, we obtain the hierarchy 
of equations as follows 



(4) 



Here ui^^^ = ii^(r, r; (?!)(n, i) ±q^) - ^^^(r, r; 0(n, t)), 



Ml =0, 



(5a) 



M- 



- u 



(5b) 



d(j) d4>dT 



(5c) 



V = d/dr, aj ~ |aj|, and Xm is the mth coordinate of the vector r, r = XmSLm/cim- 

For further consideration we have to specify the effect we are looking for and this will determine the form of 
lowest-order (j = 1) solution of Eq.(|^). Namely we will be interested in the weakly nonlinear modulation of a lattice 
wave originated by the interaction between a long wave-length acoustic mode and a high frequency mode. Thus we 
choose 



ui ^ Ao(r, t) + {^i(r,T) exp[i0(n,t)] +c.c.}, 



(6) 



where the real function Aq stands for a mean motion induced by the oscillatory wave packet, which has the complex 
envelope function Ai, and c.c. denotes corresponding complex conjugate term. Then 

u[^^^ = [exp{±iqj) - l]^ie"^("'*) + c.c, 

and in the first order [see Eqs. (^, (|5^)] we immediately arrive at the dispersion relation of the underline linear lattice 

u;2^[c^(q)]2 = 2^J2,(l-cosg,). (7) 



Next we take into account that 



duj 1 
dq uj 



[J2jSin{qj)aj/aj 



(8) 



which is the group velocity of the linear wave. Then, subject to assumption (|^) the second order equation of system 
(0) takes the form 



L 



U2 = 2*c. {((v, - v) • V) (Aie^* - ^le-'^) + x''^ (A^e^'^ - A^e'^^^)} 



(9) 
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where 



/(2) 



UJ 



-(cosgm - 1) sing„ 



(10) 



is the effective quadratic nonhnearity. 

The solvability condition for the system (in other words the conditions of the absence of secular terms in U2) 
means the orthogonality of the right hand side of Eq.(^ to the kernel of the operator L, i.e. to (H). Hence the r.h.s. 
of Eq.(^ must do not contain the terms proportional to exp(±i(/)) and we conclude that v = Vg, i.e. v introduced in 
(^ is nothing but the group velocity of the high frequency carrier wave. Next we can look for the solution U2 (it must 
be orthogonal to the first order approximation, i.e. to the kernel of the operator L) in a form of the expansion over 
the eigenfunctions of the operator L. Having done this one ensures that the only nonzero term of such an expansion 
is given by 

2^-^(2) 



Formula (|11|) is valid unless the condition uj{2c\} ~ 2(jj(q) is satisfied. As it is evident this is the condition of the 
resonant second harmonic generation . It can be satisfied in a lattice with a complex cell, but it is not difficult 

to ensure that uj{2q) 7^ '2,uj{q) for all q in a monoatomic lattice with the nearest neighbor interactions. 

Passing to the third order of the multiple scale expansion we introduce the (symmetric) group velocity dispersion 
tensor (GVDT) by the formula {v.j — duj/dqj) 

^ij = —[J^j cos{qj)aiajSij - ViVj], (12) 

UJ 

and the (symmetric) effective GVDT ilij 

^ij = — [J2jaiajSij - ViVj]. (13) 

UJ 

Then the solvability condition for the third order terms gives rise to the closed system of equations for Aq and 

l,m m 

+ ^ E a^^^ - + E (15) 

Lm m 



where 



(^m = ^^J3m(l - COSq™), (16) 
UJ 

X = — ^[2aJ3„(l - cosqrn) sing™ + 3J4m(l - cosq™)^]. (17) 



We caff Eqs.(14) and (15) the ND GDS equations. 
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III. GENERALIZED DAVEY-STEWARTSONS 



Let us now focus our attention on a special case of a 2D lattice (i.e. r = {xi,X2))- For the sake of simplicity 
the lattice will be considered symmetric, Jaj = Ja {ex. = 2,3,4 and j = 1,2) and orthogonal: ai • a2 = 0, with the 
lattice constant equal to unity, |aj| = 1. In order to diagonalize the effective GVDT Vlim in a general case we rotate 
the original Cartesian system (with the coordinate basis (1, 0) and (0, 1) ) to a new one with the coordinate basis 
ei = (Ai, A2) and 62 = (— A2, Ai), where 

"g V sm q\ + sm 172 



and Vj is the j-th component of the group velocity defined in (||) (as it is evident • = bij). In this way one 
of the directions of the new basis, namely ei coincides with the direction of the group velocity of the carrier wave, 
i. e. = VgB\ the other direction is orthogonal to it. As a result, x\ and in Eqs.(14) and (15) take the form 
x\ = e(ni — \\Vgi) and xi ~ (.(jii — \2Vgt), and envelope equations (|lj) and (^5|) are reduced to 

an^ + = -2 ( P,^, + pS) \A,\^ (19) 



drf \ dS^ drj 



i^+CA^^A, ( (3,^ +13^^] A^ + x\Ai\^ A,, (20) 



dr \ dS^ drj 

where 



5^ 52 92 

£ = 711^+722^+712^, (21) 



^ = r • ei = Aixi + A2a;2 = e(Aini + A2n2 - Vgt), 
■q = v ■ e2 = -\2Xi + A1X2 = e(-A2rii + Ain2), 

aii^—(J2-v), 022 = — , 

2J3 2J3 
/3i = [Ai(l - cosgi) + A2(l - cos(72)], lii = [Ai(l - co^q2) - A2(l - cosgi) 



1 



711 = + •^2(Af cosgi + A2Cosg2)], 722 = ^(Aa cosgi + Af cosg2), 

712 = — AiA2(cos(72 - cos 91), 



2uj 



X = — {2J3a[singi(l - cos 91) + sin 92(1 - cos 92)] + 3J4[(1 - cosgi)^ + (1 - cos 92)^]} , 

_ 4J3[sin(7i(l - cosgi) + sin 92(1 " 00352)] 
4[c.(q)]2 - [c,(2q)]2 " 



(22) 
(23) 

(24) 
(25) 
(26) 
(27) 
(28) 
(29) 



Equations ( |l9|) and (|20|) represent a generalized form of the conventional DS equations. They include the 
dispersion, diffraction and nonlinearity of the system. One of their important features is that there exists a coupling 
between the mean field (denoted by ^0) and the envelope of the carrier wave (denoted by ^1). The mean field ^0 
generates a strain field in the system. If J3 = 0, a case for a symmetric interatomic potential, we have Aq = thus the 
mean motion and hence the strain field vanishes. Another important feature for Eqs.(p^) and (|^ is their property 
of anisotropy. For different wave vector q = ((71,(72), the coefficients of the equations take different values and some 
of these coefficients may become vanishing for some particular directions of q. 

The conventional DS equations were derived firstly in surface water waves [ p+l and now are a well-known 2D 
soliton model in the soliton theory Note that for water waves, the system is isotropic (i.e. it possesses a 

continuous rotation symmetry). The envelope equations are the same for all propagating directions of the waves and 
hence the coefficients appearing in the equations are independent on qi and (72 and correspondingly /32, and 712 vanish 
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pj (see also Ref. p^). However, for the lattice system the modulating equations take a more general form because the 
lattice is anisotropic (without the continuous rotation symmetry). We mention that although the coefficients aij are 
both positive, signs of the coefficients jij may change depending on the choice of the wave- vector in the ffist Brillouin 
zone. 

We now discuss several particular cases for the 2D GDS equations derived above. In the following circum- 
stances (i.e. in some special points and lines of the Brillouin zone, see Fig. 1) the 2D GDS equations reduce to the 
conventional DS equations: 

1. qiq2 = (then A1A2 = 0, and /32 = 712 = 0) 

2. qi^q2 = q (then Ai = A2 = 2-^/2 and (32 = I11 = 0) 

3. q\ = -q2 = q (then Ai = -A2 = 2^1/2 ^^^^j ^ ^ g) 



More precisely, since an > and 0:22 > for any q, at 711722 < Eqs. I \l9\ ) and (|20|) can be classified as DSII 
equations, while for 711722 > they form a dynamical system that can be identified neither with DSI nor DSII 
equations appearing in the theory of water waves (see e.g. Jl2t). 

In the case of pure quadratic potential, J3 = 0, we have that the evolution equations for Ai and ^0 are decoupled. 
Then the GDS equations reduce to a generalized 2D NLS equation (i. e. the NLS equation plus a cross-derivative term 
Ai / (dS^drj)) . Finally, if 52 = and d/dr] = 0, the 2D GDS equations ([l9|) and (|2(]| ) recover the envelope equations 
derived in Ref. which gives rise to standard ID lattice solitons [p|,|10[. 

In the 3D case, Eqs. ( |l9| ) and (|2^) are replaced by the 3D GDS equations: 

"ll^ + "22^ + a33^ = -2(^/3l^+/^2^+/33^j l^il : (30) 

dAi , d'^Ai , d^Ai , d'^Ai / , 92 Q2 , 92 X 

I— h Ti 1 ^ „ -f Tnn ^ „ -I- Tqq —r^ -f Ti n -^r—: h Toq - — — + _ „ 



dr ' de ^^^^ 9772 +^33 Q^2 + 1^^12g^g^+ 723 5^3^+731 J ^1 

= A^ (^I3[-^ + /3^|- + /3^^) Ao + xWA,, (31) 

where a'jj, 7^'^- (j = 1, 2, 3) and x' are constants dependent on q = (gi, (?2, 93) and the parameters of the system, 
which are not needed here and not written down expliciltly. The definitions of ^, and C are given by 

^ = e(Aini + A2712 + Aana - Wgt), (32) 
7? = e(-A2ni + Ain2), (33) 
C = e(-A37ii + Xms), (34) 

where uj and A, are defined by (0) and (|l|). 



IV. LONG- WAVELENGTH LIMIT 



Note that the envelope equations (|l^ and (^5|) are invalid for q = since in this case there is a divergence in 
their coefficients. From the physical point of view this happens because vanishing q corresponds to a long wavelength 
acoustic mode in the lattice. In this case a different asymptotic expansion must be used to obtain divergence-free 
envelope equations. For simplicity we consider the case of a symmetric 2D square lattice. In this situation the 
asymptotic expansion (|^) must be replaced by 

u — uq + eui + £2^2 -I- • • • , (35) 

with 
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Uiy = U,y{^,T],T) = 0, 1, 2,- ■ 

C = e{\ini + A2n2 - ct), 
V = e^{-\2ni + A1712), 



(36) 
(37) 
(38) 
(39) 



where c = is the speed of sound, A; {I = 1, 2) are determined by the solvabihty conditions required at 0(e^)-order. 
A solvabihty condition in the forth-order of the expansion yields the generalized KP equation: 



d_ 



dv 



, d'^V J; 



dv 3J4 



, dv 



2c 



2 977^ 



= 0, 



(40) 



where v = duo/dS,. In deriving Eq.(|40|) we have assumed that J3 — eJ^ with J3 of order unity. The parameters 
A;(Z = 1,2) are direction-dependent and we find that their values can be obtained by using (|l8|) but taking the limit 
q — ^ 0. Thus the values of the coefficients in Eq.(^o|) are dependent on the ways of q approaching zero. For instances 

1. Ai - 1, A2 - 0, if 92 = 0, gi 0; 

2. Ai = A2 = I/V2, if gi = 92 = 9 ^ 0. 

The reason for appearing different values of the coefficients corresponding to different directions is also due to the 
anisotropy of the system. It is easy to see that the KP equation obtained in Ref. M is our particular case with quartic 
nonlinearity being absent (i.e. J4 = 0). Eq.(|40|) admits solitary wave solutions p^ . 

It is relevant to mention here that the coefficient of the term d^v/drf^ is positive which means that the line (i.e. 
77-independent) solitons of Eq. (BH) is stable while this equation does not admit any kind of lump (i.e. decaying when 



77- 



0) solution. 



In the same way, in 3D case Eq.(40) is generalized to 



d dv d^v 

d^ [97 ^"'9^3 



a2V 



dv 
d^ 



asv 



,dv 



a4 



duf 



as 



0, 



(41) 



where ^, 77 and ^ are the same as (|3 
by fsl) with q 0. 



(|34|). ai {I — 1, 2, 3, 4, 5) are real constants dependent on Xj {j — 1, 2, 3)(given 



V. MODULATIONAL INSTABILITY OF A PLANE LATTICE WAVE WITH A MEAN MOTION 



In recent years, the use of nonlinear envelope(or amplitude) equations for studying the stability of patterns and 
waves in systems in and outside of equilibrium is widely employed [T^JT^ . The modulational stability of a plane water 
wave (e. g. a uniform Stokes wave) was analyzed by Davey and Stewartson based on the DS equations they derived 
[ pi] . In the same way the iVD GDS equations ( p^ ) and ( p5| ) obtained here can be used to study the modulational 
stability of a uniform Stokes lattice wave in N dimensions. A Stokes lattice wave here means a linear plane lattice 
wave with the wave vector q. 



(42) 



Note that the uniform vibrating solution of Eqs.(|14|) and (|15|) reads 

Ao = 0, = [/oCxp(-ifiT), 



which, when incoperating the carrier wave (see (^ ), corresponds a plane lattice wave with the wave vector q and the 
frequency w(q) + excited in the system, where Uq is a constant and ft = xUq- Assume that a perturbation is added 
into the uniform vibreting solution (|42|), i. e. 



Ao{xi, X2, t) = K+exp(i^(3, 

m 

Ai{xi, X2, t) = ?7o exp(-ifJT) 



m 

1+6+ exp(i ^ QmXm) + £- exp(-i ^ Q,, 



(43) 
(44) 



7 



with k±{t) — K±(0) e'xjp[{aji ± i<Ji)T] and £±(t) = e±(0) exp[(crjj ± Z(T/)r], where Q=(Qi, Q2, Qn) and (T/=cr/(Q) 
are respectively the wave vector and frequency of the perturbation, an = (Tfl(Q) denotes the growth rate of the 
perturbation, k±(0) and £±(0) are small constants with the condition k1(0) = k+(0) because of the reality of ^o- 
Substituting (Esl) and into the Eqs.dlJ) and (|l5|) we obtain a set of linear equations on k±(0) and e±(0): 



- (aiiQ? + a22Ql)K+{0) + 2zt/2(/3iQi + /32O2)[£+(0) + £l(0)] = 0, 
{aiiQl + a22Ql)n-{0) + 2^C/2(/3lQl + l32Q2)[e-{0) + e*+m = 0, 

{n + ia- 7ng? - J22QI - 712Q1Q2 - 2xC/o']£+(0) - xC^o£l(0) ^ *(/3iQi + /32Q2)«;+(0) = 0, 
(r! + la* - -fnQl - I22QI - 712Q1Q2 - 2xUl]e-{0) - xt/o^W + ^{PiQi + /32Q2)'*-(0) = 0, 



(45) 
(46) 
(47) 
(48) 



where a — aR + iaj. A solvability condition of Eqs.([45|-[48D results in 




-x + 



El. 

rr 



(49) 



Note that the right side of Eq. (M) is real. Thus when 




%QlQr. 




2 (Em '^mQm)' 

r. 



4 ^ 

l.m 



^ImQlQm / > 0, 



(50) 



one has ct/ = 0. As a result if the condition 150) is satisfied we have the growth rate 





2 (Em ^" 



E, QiQ n 



1/2 



(51) 



Thus one always has a positive an branch if the condition (|5C|) is satisfied. In this case the perturbation grows 
exponentially and hence the unifrom vibrating solution is modulationally unstable. 



For the 2D GDS equations (|19|) and (|0|) , the condition of the modulational instability (|50| ) reads 

2{PlQl+P2Q2?' 



{111QI + 122QI + 112QIQ2) 



"22(51 



- ^(7ii< 



722*52 +7l2(3l(32) \ > 0. 



(52) 



Thus due to the anisotropy of the lattice (i. e. ^2^12 7^ 0) the criterion (^2|) gives much richer behavior for the stability 
of the Stokes wave than that in isotropic systems (e. g. water waves). In particular, for a given Stokes lattice wave 
there exist two (or may be four depending on the Stokes lattice wave) wave vectors Q for which the instability evolves 
with the biggest incriment. This phenomenon reminds the so-called strengthing of inhomogeneities, known in the 
theory of beam propagaion in Kerr medium p9[ | . There is however an essential difference originated by the anisotropy: 
the biggest exponent is characterized by the amplitude of the value of the wave vector and also by the lattice direction. 
The position of the points providing the largest incriment depends on the choice of the wave vector of the Stokes 
lattice wave. 

The outcome of this type of instabilty may result in the fromation of solitons or the appearence of homoclinic 
structures(see Sec. 3.3 of Ref. |12|). 



VI. SOLITON SOLUTIONS 



We now consider the soliton solutions of the nonlinear evolution equations derived above. Taking 2D GDS equa- 
tions (p^ and (|2^) as an example, to obtain the soliton solutions we employ Hirota's bilinear transformation method, 
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an ingenious technique of finding exact multi-soliton solitons for nonlinear evolution equations ]17| , |18| . Introducing 
the dependent variable transformation 

^o = -4(^/3i^+/?2|^)logF, A^^G/F (53) 

with F (real) and G (complex) being the functions of t, ^ and 77, Eqs.(p^ and ( |20| ) arc transformed into the following 
bilinear form 

{aiiDl + a22Dl)F-F=\G\'', (54) 

{iDr + 711^1 + + li2D^D,,)G • F - 0, (55) 

[(711 - 2Pl)Dl + (722 - 2/3|)A' ^ ^(^^^ _ 4/3i/32)i?5A,]F • F + x\G? - 0, (56) 



where D^, -D^ and are Hirota's bilinear operators defined by |[l7|Jlq | 

G{^,V,r)F{e,v',r%,=^,^,=„,r'=r. (57) 



DJ'Dl^DPG ■ F 



t)777 lySr 9r' 

In order to get a one-soliton solution we assume 

F = l + Lcxp($ + $*), G = exp($) (58) 

with 

* = (Pi? + ipi)^ + (to + iqi)v + {sB. + isi)T + $ofl + «*o/, (59) 



where L,pif ,pi,q]i,qi, sh^ si.^on and <I>o/ are real, yet to be determined constants. Substituting (58) into (p34|)-(|56|) 
we obtain the set of algebraic equations 

8L{aiipl + a22ql) -1 = 0, (60) 

711 b/?, - v]) + luipRm " piqi) + 722 ((jI - qj) - s/ = 0, (6I) 

2(7iiPflP/ + l22qRqi) + li2{pRqi +PiqFi) + sr = 0, (62) 
X + 8L[(7ii - 2pf)pl + (722 - 2/3|)g^ + (712 - 4/3i/32)Pflte] = 0. (63) 



L = 77, 2 2T- (64) 

8(aiiPfl + "22'?^) 



From Eq. (pOD we get 



Eqs.(61) and (|62|) give rise to the "dispersion relations" 

Si = iiiip^i- p]) + ii2ipRqR - Piqi) + 722{q% - qj), (65) 

sr = -2(7iipflp/ + 7229^9/) - iMPRqi + PiqR), (66) 



with pr^pi, qn and q/ being arbitrary constants. Eqs.(63) gives a condition for the one-soliton solution. 
From ( 53 ) and the results given above we have 

Ao = -A{f3ipR + /32te)[l + tanh(0 - So)], (67) 
Ai ^ [2ianpl + a22ql)]^/^sech{9 - So) exp(i(p), (68) 

with 9 = PB.S. + qRV + srt + $ok, f ^ Pi£. + qiV + sjt + $0/ and 6q = (1/2) log[8(aiip|j + a22g|)] (^'oi? and $0/ 
are arbitrary constants). Thus the single-soliton solution obtained is a line soliton, which consists of two parts, a 
vibrating wave packet (Ai, a envelope soliton) and a mean displacement field (Aq, a kink). 
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The two-soHton solutions of the Eqs.(|T9|) and ( pO| ) can be obtained by chosing 

F = 1 + Li exp($i + <Pl) + L2 exp($2 + ^'2) 

+ (L3 + iLi) exp($i + + (L3 - iLi) exp($J + $2) + L5 exp($i + $2 + + ^-2), (69) 
G = cxp($i) + exp($2) + (Afi + iM2) exp($i + $2 + *I) + {M3 + iM^) exp($i + $2 + ^l), (70) 

with $j = (pji?+«Pj/)C + ('7jfi + »'?j7)7?+(sjfl + iSj/)T + $°j^ + i$°^ (j 1, 2), where p^-fl, p^j, qjR, q^j, SjR, Sjj, 
and ^'jj are real constants. When (^) and (^0|) are substituted into the bilinear equations (|5^)-(|56|) we obtain a set 
of nonlinear algebraic equations for the real coefficients Lj {j = 1, 2, 5) and Mj {j = 1, 2, 3, 4) appearing in ( |69| ) 
and (^. Solving these equations one can get the expressions of Lj and Mj, as well as the "dispersion relations", 
SjR,i — SjRj{pjR,Pji , qjR, qji) {j — 1, 2), which have been given in the Appendix A. To guarantee ( |69|) and (|70| ) are 
two-soliton solution, the following conditions must be imposed 

712 = 4/3i/32, (71) 
Qii ^ «22 ^ _ ]_ ,„2^ 

711-2/3? 722-2/32 ^ I ; 

In addition, for p2B. and 52/?,, there is a constraint 

a22(aiiP2fl + a22g2fl)x = ("ii/?! + a22/3?)P2fl + 2a22/3|9L. (73) 

It is easy to show that the integrable conditions of the standard DS equations (i. e. the ones amenable to be solved 
by the inverse scattering transform) derived in water wave problem are the particular case of the conditions ( [zif ) and 
(i'2) (see Appendix B). This fact implies that the GDS equations (^ and ( pO| ) may be integrable under the conditions 
© and @). 

We note that different equalities in these conditions, however, reflect different physical properties. In particular, 
J) and the first equality in ( [7^ ) result in an equation for the wavevector only (i.e. having the form / (91,92) = 
where f{qi, (72) does not depend on the lattice parameters, i.e. on J2) for which the existence of solitons is possible. 



Then first two equalities in (71) and ( [72| ) allows one to find the particular values of the nonlinear coefficients. In other 
words, the above conditions specify the set of points in the first Brilloun zone, and necessary values of the nonlinear 
forces. What is important for the next consideration, that such points in the Brillouin zone do exist. Indeed, as an 
example we mention that the above conditions are satisfied for all points q = (gi, 0) and q = (0, 52)- 



Eqs. ( |69| ) and (70) describe two obliquely interacting solitons in the (^,?y) space. The interaction results in a 
phase shift (i. e. position shift) for each soliton. 

It is possible to get iV-soliton solutions of the 2D GDS equations (|l^) and ( pO| ) using their bilinear representation, 
Eqs.(54)-(p6|) under the integrable conditions ( [7l| ) and (|7^. We note that due to the anisotropy inherent in the lattice 



system (i. e. 7^ 0), the existence of the two-soliton solution needs the condition 712 = 4/3i/32 (Eq.(|7l|) ), which is 

absent for isotropic systems (e. g. water waves). 



VII. CONCLUSION 



Using a quasi-discrete multiple-scale method we have derived the envelope equations of weakly nonlinear modula- 
tions of A^-dimensional lattice waves. The equations are obtained for the case of interaction of a highly frequency mode 
with a long-wave-length acoustic one (also called mean field) and can be classified as generalized Davey-Stewartson 
equations. In the case at hand, due to the anisotropy of the lattice system, the GDS equations in two dimensions are 
reduced either to the DS equations or to a form which does not appears in the theory of water waves [Til ]. The mean 
field coupled to the oscillatory short wave packet results from the cubic interatomic potential in the lattice. Addi- 
tionally, generalized Kadomtsev-Petviashvili equations describing the evolution of a long wavelength acoustic mode 
in the lattice are also presented. We have also studied the modulation instability of Stokes waves and provided some 
exact soliton solutions for the two-dimensional GDS equations based on Hirota's bilinear transformation method. 

The results reported here recover the known ones in one-dimensional systems, which give rise to standard lattice 
solitons. On the other hand the method can also be used to study the weakly nonlinear modulations of the wave 
packets in vector lattices or in lattices with a complex cell. The derivation procedure involves more cumbersome 
calculation but the envelope equations obtained take still a form similar to ( |l^ ) and (|l5| ) for high frequency wave 
packets and (^) and (^ for long wave acoustic modes. 
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Appendix A 

The expressions of Lj and Alj for two-sohton solutions appearing in ( |69| ) and ( |70| ) are given by 

8(aiip2^ + a22qlii) ' 

= I , 

8(aiiP2_R + "2291^) ' 

, _ 1 aiiri+ + Q!22EI^_ 



2 [aiir:+ + a22Sl+]2 + 4[anA^+ + a22Al+]2 ^ 
ailA'^_l_ + q;22A'L^ 
^ [anri+ + a22El+]2 + 4[anAP + + a22A%]^ ' 

1 An 



64 A 



5d 



8 Mid' 2M2d' 8 Mad' 2 M4d ' 

L5n = a?i(rt_)2 + a22(St_)2 + 2aiia22(r:_E:_ + 4A^_Al_), 

L5d = (aiiP?_R + a22g?i?)(aiiP2fi + "22g|_R,) • 

[a?i(ri+)2 + aU^t+r + 2ana22(El+ri+ - 4AP_+A1+)], 

+c^22{i'liR - <l2Rf + (91/ - q2if[{qii - q2if - '^{Sqfjj - 92fl)]} 
+2aiia22{[(pi/ -P2/)^ - + pD] [(91/ - <l2i f - {qln + qln)] 

-A{pii - P2l){qu - q2l)(piRqiR. - P2R.q2B) - 'iqiB.q2RPlRP2R} , 

Mid = M2d = + a22qlR)[ali{Tt+f + a22(Et+)' + 2ana22(r:+E:+ + 4AP +A1+)], 

M2n = "nPlflbl/ -P27)[(P17 -P2/)^ ^ PlR + pIr] + a22«li? (i?!/ - 92/) [(^l/ " 92/)^ " 'Z?/? + 927?] 

+aiia22{(gi/ - g2/)[9ii?ri_,_ + 2piRP2R{qii + q2R)] + [pii - P2i)[piR.^Z^ + 2qijiq2RP2R]), 
M3n = " A,)' + {Pii - P2if[ipii - P2if - 2{3pIr-pIi^)]} 

+a22{(9iK - 92fl)^ + (91/ - 92/)^[(9i/ - 92/)^ - "^iSqln - Qir)]} 
+2aiia22{[(pi/ -P2/)^ - (pIr + plRMiqii - q2if - {qln + 9^)] 

-4(pi/ - P2l){qil - q2l){PlRqiR - P2_R92i?,) - 'iqiRq2RPlRP2R}, 

Msd = Mid = HipIr + a2292fl)[«?i(rt+)' + "22(2!+)' + 2ana22(r:+EI+ + 4A%Al+)], 

M4„ = -aliipii - P2i)P2r[{pii - P2if +pIr-pIr] ~ a22(9i/ - 92/)92i? [(91/ - 92/)^ + 9i7? - 92ii] 
-aiia22{(<7i/ - 92/)[92/?ri_|_ + 2pijip2R{qii + q2R)] + [pii - P2i)[P2R^Z+ + 2qiRq2RPiR\}, 



where 



^tia2 = (PlI + ^lP2lf ± [PlR + Cr2^'2i^)^ 
^^i<T2 = (91/ + Cri92/)^ ± {qiR + CT292i?,)^ 
^lia2 = (PlI + ^lP2l){piR + Cr2P2R), 
^lia2 = i^lil + 0"l92/)(9li?. + 0-292i?), 



with 0-j = ±1 {j = 1,2). 

The "dispersion relations" are given by 



Sir 



SlI 



-4/3i/32(pifl9i/ + PiiqiR) + 2/3| ( '^"^^^^^'^ - qiiqiR 

V "22 



1 

"110:22 



-{[-"?i/3|b?fl - pIi) " «22/3f (9?fl - 9u) 



2/32 



/ "2291/91/? 



"11 



Pl/Pl/? 
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+Q!iiQ!22[/?i(Pifl -p?/) + 4/?i/?2(piflg'ifl -Pi/gi/) + - 9?/)]}, 

S2R = -4^PlP2{P2Rq2I + P2iq2R) + ^P2 q2iq2R +2/3^ P2IP2R 

V "22 / V an 

= r. \ {[-^11/^2 (P2fl - pIi) - <A21^1{<iIr - ill) 

+Q!llQ!22[/?i (P2ii - P\l) + ^0102iP2Rq2R - P2iq2l) + P'HqlR - qh)]} ^ 

where pjn, pji, qjn and qji [j = 1, 2) are arbitrary constants. 
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Appendix B 

One type of the standard DS equations which can be solved by the inverse scattering transform is (see p. 240 in 
Ref.[12] for the case of r = —q*) 

^ ^ = ~2— (\q'\) 

dx^ dy^ dx^ ' 

.dq 1 2 cl'^l 1 9^1 , I i2 

with cr^ = ±1. Taking the transformation x y rj^ t 3^'''^ 1 ~* "71^^' ^ ^ above equations 

become 

Comparing with Eqs. (|l^) and (^0|), for the last two equations we have 

an = 0-2, a22 = -1, /9i = -1, /32 = 0, 
711 = 1. 712 = 0, 722 = cr"2, 

which satisfy the integrable conditions (fzi] ) and ([72|). 
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